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Abstract 

In this review, based mainly on [3], we investigate a key information-theoretic per- 
formance metric in multiple- antenna wireless communications, the so-called outage 
probability. This quantity may be expressed in terms of a moment generating function, 
for which we present two separate integral representations, one involving a particular 
a— form of Painleve V. We also present a representation involving a non-linear second 
order difference equation. 

1 Introduction 

In Random Matrix Theory (RMT) many statistical quantities can be described as deter- 
minants. This is especially true in the cases of what are known as the classical ensembles, 
for example, the Gaussian Unitary Ensemble (GUE) or the Laguerre Ensembles. Since the 
groundbreaking work of Tracy and Widom [20], which characterized the largest eigenvalue 
distribution for the GUE, one common approach is to write some statistical quantity as a de- 
terminant and then express the determinant as something involving a Painleve transcendent, 
a solution to one of the classical Painleve second order non-linear differential equations. 

Techniques to find the determinant and then the resulting differential equation are quite 
complicated. Tracy and Widom used more of an operator theory approach, others have used 
an integrable system approach, whilst others still have used a stochastic equation approach. 

This paper highlights a technique known as the "ladder operator" approach and is, in 
some sense, one long example illustrating the technique and describing the application of 
interest. The specific application relates to the performance of multiple- input multiple- 
output (MIMO) wireless communication systems, in which both the transmitter and receiver 
devices are equipped with multiple antennas. Such systems have been the subject of intense 
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interest since the key papers [T7] and [5], and now form the cornerstone of most modern 
day wireless systems (Wi-Fi, cellular networks, etc). Here, a key fundamental performance 
measure is considered, the so-called "outage probability", and this is shown to involve the 
probability distribution of a certain "linear statistic" in a Laguerre random matrix ensemble. 
The problem thus falls naturally within the realm of the ladder operator framework. 

Here is the idea of the approach. For the example at hand, the linear statistic of interest 
can be described through its moment generating function as a Hankel determinant, which 
using the theory of orthogonal polynomials (for a certain nonstandard weight generally) 
can be computed via a product of norms of monic orthogonal polynomials. Now it is well- 
known that orthogonal polynomials satisfy three term recurrence equations. In fact, the 
two coefficients in the recurrence combined with initial conditions, completely determine the 
polynomials. Thus information about the coefficients in the recurrence equations should 
yield information about the Hankel determinants. 

The path to this information is from the ladder operators, two formulas that connect the 
polynomials one index apart to their derivatives. These yield, using basic complex analysis, 
a set of equations in the coefficients along with two additional auxiliary quantities that arise. 
The story would end here, except for the fact that often there is a "time" parameter implicit 
in the original weight and thus in the polynomials themselves. Using "time" evolution one 
can then, using only elementary means, find a pair of coupled Ricatti equations in the two 
auxiliary quantities. These then lead directly to a Painleve equation. It should be pointed 
out that this method works at least in principle, if the "time" parameter is present and if 
the derivative of the logarithm of the weight is a rational function. Then one can in many 
cases follow the steps illustrated in this paper. 

Here is an outline of the paper. The next section contains the preliminaries of the theory 
including the ladder operator equations. Section III shows how the application of interest, 
the outage probability performance measure which arises in the application of MIMO wireless 
communication, can be described using the RMT framework. Section IV provides the details 
of the path to the differential equation solutions, which are presented Theorems 1 and 2. 

2 Preliminaries 

2.1 Linear Statistics of Hermitian Random Matrices and Hankel 
Determinants 

For the MIMO capacity application, it will be seen that the problem of interest falls within 
the general theory of linear statistics of Hermitian random matrices, with a close connection 
to the theory of orthogonal polynomials. Here a brief introduction to the general theory is 
given, and preliminaries are established for later use. 
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We will require the distribution of a certain linear statistic 



N 



(2.1) 



k=l 



in the eigenvalues {x k } of a N x N Hermitian random matrix, with joint eigenvalue density 
of the form 



N 



p(x 1 ,...,x N ) (xY[w {x k ) Yl ( x j- x i) 2 > x k e (a,b) 

l<i<j<N 



(2.2) 



k=l 



for some weight function ioq(-) . It is convenient to attempt to characterize the distribution 
of the linear statistic (12. ip through its moment generating functional], 



M(X) =E 



N 



exp A ^2 f( x k) 



k=l 



E 



' N 

n 

,k=l 



(2.3) 



which upon substituting for (12.21) gives 

M f(a,b) N ril<i<j<A l ' ! I 



M{\) 



■>U - x i) 2 rifc=i w(x k )dx k 



N\ f(a,b) N ril<j<j<A r (- E i X *) 2 11^=1 w o( x k)dXk 



where 



w(x) := w (x)e Xf( - x ^ 



(2.4) 



denotes the deformed version of the reference weight wq(-) . Application of the Andreief- 
Heine identity [16] now directly leads to 



M(X) 



D N [w] 



det^/^x 1 " 1 "- 7 2 w{x)dx 



N 



Dn[w ° ] det (j b a x^- 2 w (x)dx) N 

V a / 1,7 = 1 



(2.5) 



which is a ratio of Hankel determinants. By virtue of the Selberg integral, for most "classical" 
weight functions of interest, the Hankel determinant in the denominator of (12 .4p admits an 
explicit closed-form (non determinantal) representation. The numerator, on the other hand, 
is much more difficult to characterize, since it involves the more complicated deformed weight 
w . To proceed, methods based on orthogonal polynomials will be introduced in the sequel. 



lr rhe parameter A is an indeterminate which generates the random variable 53fe=i f( x k) 
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2.2 Orthogonal Polynomials and their Ladder Operators 

We start by noting that 

II - = det ( ''} ^ ; = det (Pi-l(^))Z=l (2-6) 
l<i<j<N 

where Pj(-) represents any monic polynomial of degree j , 

P 3 (z) = zJ + Pl (j)zJ- 1 + ... (2.7) 

Applying this in (12 .4p and once again integrating via the Andreief- Heine identity, the nu- 
merator evaluates to 

D N [w] = det ( I P i ^ l (x)P j ^ 1 (x)w(x)dx\ . (2.8) 

\Ja J i,j=l 

If we orthogonalize the polynomial sequence {P n (x)} with respect to w(x) over the interval 
[a, b], i.e., 

b 

Pj(x)P k (x)w(x)dx = hj6 j>k , j, k = 0, 1, 2, ... (2.9) 
with hj denoting the square of the L 2 norm of Pj over [a,b], then (12.81) reduces to 

N-l 

Dn[w] = h k . (2.10) 

k=0 

The key challenge is how to characterize the class of polynomials which obey the orthog- 
onality constraints in (12. 9p or, more importantly, the norms of such polynomials required to 
evaluate f[2~TU|) . 

If all the moments of the weight w exist, then the theory of orthogonal polynomials 
states that the P n (z) for n — 0, 1, 2, ... satisfy the three term recurrence relations, 

zP n {z) = P n+1 (z) + a n P n {z) + n P n - X {z). (2.11) 

The above sequence of polynomials can be generated from the orthogonality conditions (12.91) . 
the recurrence relations (12.111) . and the initial conditions, 

P„(*) = l, PoP-i(z)=0. (2.12) 

Substituting (12. 7p into the recurrence relations, an easy computation shows that 

p x (n) - p 1 (n + l) = ain, (2-13) 
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with p 1 (0) := 0. A telescopic sum of (12 . 1 3[) gives 



n-1 



3=0 

From the recurrence relation (12.111) and the orthogonality relations (12. 9p . we find 

Pn = T~~~ ' (2-15) 

n n-l 

We shall see that Pi(n) plays an important role in later developments. For more infor- 
mation on orthogonal polynomials, we refer the reader to Szego's treatise [15] . 

Next, we present three Lemmas which are concerned with the "ladder operators" asso- 
ciated with orthogonal polynomials, as well as certain supplementary conditions. Note that 
these have been known for quite sometime; we reproduce them here for the convenience of the 
reader using the notation of [2J, where one can also find a list of references to the literature. 
We also mention that Magnus [9] was perhaps the first to apply these lemmas — albeit in a 
slightly different form — to random matrix theory and the derivation of Painleve equations. 
Tracy and Widom also made use of the compatibility conditions in their systematic study 
of finite n matrix models [19] . See also [3]. 

Lemma 1 Suppose v(x) = — \ogw(x) has a derivative in some Lipshitz class with positive 
exponent. The lowering and raising operators satisfy the differential- difference formulas: 

P^z) = -B n (z)P n (z) + (3 n A n (z)P n ^(z) (2.16) 

P' n _ x {z) = [B n (z)+y'(z)}P n -x(z)-A n ^(z)P n (z), (2.17) 

where 

A n (z) := J_ I A*)-AV) P *(y) W (y) d y ( 2 .1 8 ) 

K J a z-y 

I f b y'(z) -v'(v) 

B n {z) := / -U ^-P n (y)P n _ 1 (y)w{y)dy. (2.19) 

K-i J a z ~y 

A direct computation produces two fundamental supplementary or compatibility conditions 
valid for all z G C U {oo} . These are stated in the next Lemma. 
Lemma 2 The functions A n (z) and B n (z) satisfy the supplementary conditions: 

B n +i(z) + B n (z) = {z- a n )A n {z) - v'(z) (Si) 

1 + (z - a n )[B n+ i(z) - B n (z)] = (3 n+1 A n+ i - /3 n A n _i(z). (S 2 ) 



It turns out that there is an equation which gives better insight into the coefficients a n and 
(3 n , if (Si) and (S 2 ) are suitably combined to produce a "sum rule" on A n (z). We state 
this in the next lemma. The sum rule, we shall see later, provides important information 
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about the logarithmic derivative of the Hankel determinant. 
Lemma 3 The functions A n (z), B n (z) , and the sum 

n-1 

EM*), 

3=0 

satisfy the condition: 

n-1 

B$(z) + v'{z)B n (z) + M z ) = Pn A n {z) A n ^(z). (S' 2 ) 

3=0 



3 Information Theory of MIMO Wireless Systems 

In this section we introduce the wireless communication problem of interest, and connect it 
with the general linear statistics framework introduced previously. 

We consider a MIMO communication system in which a transmitter equipped with n t 
antennas communicates with a receiver equipped with n r antennas. Denoting the transmit- 
ted signal vector as x G C n * and the received signal vector as y G C™ r , under a certain 
assumption on the channel (known as "flat fading"), these signals are related via the linear 
model 

y = Hx + n, (3.1) 

where n G C nr the receiver noise vector, is complex Gaussian with zero mean and co- 
variance E(nxv) = I n ,. . The matrix channel matrix, H G C n ' xrat , represents the wireless 
fading coefficients between each transmit and receive antenna. The channel is modeled 
stochastically, with distribution depending on the specific wireless environment. Under the 
realistic assumption that there are sufficient scatterers surrounding the transmit and receive 
terminals, the channel matrix H is well modeled by a complex Gaussian distribution with 
independent and identically distributed (i.i.d.) elements having zero mean and unit variance. 
This matrix is assumed to be known at the receiver^], but the transmitter only has access to 
its distribution. The transmitted signal x is designed to meet a power constraint: 

£(x f x) < P . (3.2) 

Our objective is to study the fundamental capacity limits of a MIMO communication 
system. Such limits are described by the field of information theory, founded by Claude 
Shannon in 1948 [TJ]. Specifically, information-theoretic measures allow one to precisely 
determine the highest data rate that can be communicated with negligible errors by any 
transmission scheme. Consequently, information theory offers a benchmark for the design 

2 In practice, this information can be obtained using standard estimation techniques. 
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of practical transmission technologies, and has become an indispensable tool for modern 
communication system design. 

The capacity of a communication link is determined by the so-called "mutual informa- 
tion" between the input and output signals. For the MIMO model (13 .ip it is given by: 

J(x;y|H)='W(y|H)-7£(n) (3.3) 

with H(y|H) denoting the conditional entropy of y , 

H(y\U) = - [ p(y|H)logp(y|H)dy, (3.4) 

where p(y|H) denotes the conditional density of y given H. This formula represents the 
maximum amount of information that can be reliably transported between the transmit- 
ter and receiver (i.e., it represents the rate which is "supportable" by a given realization 
of the MIMO channel). It was proved in [T7] that the conditional mutual information 
7(x;y|H) is maximized by choosing the input signal vector x according to a zero-mean 
circularly-symmetric complex Gaussian distribution with covariance Q x = E (xxH satisfy- 
ing tr (Q x ) < P . In this case, the mutual information (13. 3 p was shown to be 

7(x; y|H) = log det (v + HQ X H^ . (3.5) 



In this paper, we will consider a scenario in which the channel is selected randomly at 
the beginning of a transmission, and remains fixed during the transmission. In this scenario, 
it is impossible to guarantee that the communication will be completely reliable, since no 
matter what transmission rate R we choose (which is assumed fixed) there is always a 
non-zero probability that the rate may not be supportable by the channel. In other words, 
there is always a chance that the mutual information 7(x;y|H) falls below 7?, and thus 
communicating at rate R becomes impossible. This is referred to as an "outage event" , and 
the probability of this occurring is called the outage probability, 

P out (R) :=Pr(7(x;y|H) < R) . (3.6) 



Here, we will make the common assumption that 

Qx = "In, , (3.7) 

n t 

corresponding to sending independent complex Gaussian signals from each transmit antenna, 
each with power P/n t . Hence, with this input signal covariance, the quantity P (> 0) 
will also represent the signal-to-noise ratio (SNR). With Q x given by (13.71) . the mutual 
information 7(x;y|H) becomes 

7(x; y |H) = log det f I Ur + ilfflA , t:=%. (3.8) 
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To fix notation, let M := max{n r ,n t }, N := min{n r , n t }, a := M — N and define 

_ f HH^ n r < n t 
\ HtH, n r > n t 

with positive eigenvalues {xk} k=1 • With these definitions, and with det(I + AB) = det(I + 
BA) , we can further evaluate 

N ( 1 \ 
7(x;y|H)=^log [\ + -x k J 

N 

= -Nlogt + Y,^g{t + x k ). (3.9) 

k=l 

Computation of the outage probability (13. 6p requires the probability distribution of 
J(x;y|H). From (I3.9p . this is clearly a linear statistic in the eigenvalues of the Hermi- 
tian random matrix W (with a constant shift of — iVlogt). Moreover, W is complex 
Wishart distributed [IB], thus the eigenvalues {a^}^ are well-known to admit the joint 
density 

N 

p(xi,x 2 ,...,xn) cx JJitfLagOr/) Yl ( x k~Xj) 2 , x\ 6 [0, oo) (3.10) 

1=1 l<j<k<N 

where 

w Lag (x) := x a e~ x , a > -1, x € [0, oo) (3.11) 

is the classical Laguerre weight. 

Our aim will be to compute the moment generating function of the linear statistic, 

M(X) := E u (>'( x ;yl H )) = t~ NX M(X) (3.12) 
where A4(X) is identified by (12.31) but with the following particularizations: 

(f(x),w {x),w(x),a,b) (\og(t + x),w Lag (x),w dLag (x,t),0,oo) (3.13) 
where w^^x, t) is a deformed Laguerre weight, 

w dLag (x, t) := (x + t) x w Lag (x), t>0 x>0. (3.14) 
Thus, (12. 5p immediately gives 

= (3.15) 

where 

D N {t, A) = det (/i i+j _ 2 (t, \))f ij=1 (3.16) 
8 



is the Hankel determinant generated from Wdhag{x,t) with moments 

Hk(t,\):= x k w dLag (x)dx , k = 0, 1,2, (3.17) 

■> 

The quantity D^lw^g] in the denominator of (I3.15P is the Hankel determinant generated 
from the classical Laguerre weight, w^^x) , and can be computed in terms of the Barnes 
G -function as 

„ r . G(N + 1)G(N + a + 1) . . 

£>ivKag] = -^ J {a \ 1} L , G(l) = l. (3.18) 

Our next objective will be to compute a non-determinantal representation for the (scaled) 
moment generating function (13.151) . This, in turn, will require evaluation of the Hankel 
determinant D N (t,\) in (I3.16p . We will address this problem in the sequel by appealing to 
the orthogonal polynomial framework introduced in Section 12.21 We should like to mention 
here that, unlike the classical ladder operators, the "coefficients" in our ladder operators are 
V dependent, as we shall see later. 



4 Painleve V and the Continuous and Discrete a -Form 
4.1 Main Results 

The following two theorems present the main results of the paper: 

Theorem 1 The Hankel determinant Djv(t, A) admits the following integral represen- 
tation: 

D N (t,\) NX 



/)v I t ™P{J t f(y(s),y'(s),s)dsj (4.1) 

where 

fjyjs), y>(s),s) := ^ + 2 (* + " ' % ± ^L± (f +f Z ^ ± 2 ^ ± *V (42) 

-2(2iV g + g(s + a - A))?/ 3 + a 2 y A - [g ?/( g )] 2 

4y(y-l)a ' l ' J 

Theorem 2 Equivalently, the Hankel determinant D^{t,X) also admits the following in- 
tegral representation: 

W>-«"W/' g »M-*V) (4.4) 



where H^{t) satisfies the Jimbo-Miwa-Okamoto a— form of Painleve V: 

{tH'iff = [tH' N -H N + H' N (2N + a + A) + NX] 2 

- A(tH' N -H N + S N ) [(H' N ) 2 + XH' N ] (4.5) 
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with 5 N := N(N + a + A) . 

Before presenting the proof of these results, we would like to point out that Painleve 
equations first appeared in the early 1900's through the work of Painleve and his collabora- 
tors [6]. In the mid 1970's, Painleve equations first appeared in characterizing the correlation 
function of an Ising model through the pioneering work of Barouch, McCoy, Tracy and Wu, 
see [10]. The 1-particle reduced density matrix was shown in 1980 to satisfy a particular 
Painleve V, see [TJ. For a recent review on this and other related problems in matrix ensem- 
bles, see [32]. Another Painleve V appeared in the Hankel determinant associated with the 
"time evolved" Jacobi polynomials, see [JJ. 



4.2 Proof of Theorems 1 and 2 

4.2.1 Compatibility Conditions, Recurrence Coefficients and Discrete Equa- 
tions 

For the purpose of applying the ladder operator framework introduced in Lemmas 1-3, an 
easy computation shows that 

v(z,t) := - log w dLag (z,t) = -alogz - Alog(z + i) + z, 



and therefore 



i, \ at A 

v'CM = — + l 

z z + t 



v'{z,t)-v>(y,t) a A 



z-y zy (z + t){y + t) 

Substituting the above into ( 12.181) and (12.191) . followed by integration by parts, we obtain 

A n (z)= l -^l + ^l (4.6) 

B n (z) = - n + r " (t) + r -«® (4.7) 
z z + t 

where we have introduced the auxiliary quantities: 



A r°° \P (r)\ 2 

Rn{t) := t- / [ -^-w dUe (x,t)dx (4. 

"n Jo X + t 



A f°° P n (x)P n ^(x) 



r n (t) := / — w d hag(x,t)dx. (4.9) 

"n-l JO X + t 

The next Lemma gives a representation of the recurrence coefficients a n , (3 n , Ylj anc i 
Pi(n) — the coefficient of z n ~ x of P n (z) — in terms of the auxiliary variables r n and R n . 
Note that Pi(n) also depends on t but we do not display this if there is no confusion. 
Lemma 4 The recurrence coefficients a n and (3 n are expressed in terms of the auxiliary 
quantities r n and R n as: 

a n = 2n + l + a + X-tRn (4.10) 



1 — Rn 



r n (2n + a + A) + Tn A? " + n(n + a) 



(4.11) 
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Furthermore, 



n-l 



tYj R J = n(n + a + X) - /3 n -tr n , (4.12) 

j=0 

Pl (n) = -P n -tr n . (4.13) 

Proof We start from (Si). Equating the coefficients of z' 1 and [z + 1)^ 1 , we obtain the 
following difference equations relating a n to r n and R n : 

-(2ra + 1 + r n+1 + r n ) = a- a n (\ - R n ) (4.14) 

r n +i +r n = X — R n (t + a n ). (4.15) 

To proceed further, we take note of (I4.6p -( 14T91) . and derive identities based on the supple- 
mentary condition (S' 2 ) , which will be of particular interest. A straightforward (but long) 
computation shows that the r.h.s. of (S' 2 ) becomes 

n-l 

B 2 n (z)+y'(z)B n (z) + Y J A 3 ^) 

= z~ 2 [(n + r n ) 2 + a(n + r n )j 

n-l 



+ z 1 <n 



Rj + r n [X - a - t - 2(n + r n )]/t + (n - A)/i} 

j=0 
n-l 

+ (z + Rj + r n [t + a-X + 2(n + r n )]/t - nA/i} 



i=o 

+ + i)~ 2 [^-Ar n ]. 

Now focusing on (S 2 ) as presented above and equating the coefficients of z~ 2 , (z+i) -1 , 

[z + t)~ 2 , give rise to the following difference equations involving (3 n , r n , i? n and ^ ■ Rj : 

(n + r n ) 2 + a(n + r n ) = (3 n (l - R n )(l - R n .{) (4.16) 



n - J2 Rj + j [A - a - t - 2(n + r n )] + = T [(1 " + ^ ~ ^-iW ( 4 - 17 ) 

j=0 



Rj + j[t + ol - A + 2(n + r n )] - — = [(1 - i^i^-i + (1 - fl^i)/^] (4.18) 

3=0 



r 2 n -Xr n = p n R n R n ^. (4.19) 
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Observe that (I4.17P and (I4.18P are equivalent. We shall see later that (I4.17p . when combined 



with certain identities, performs the sum Xyj=0 automatically in closed-form. This sum 
will provide an important link between the logarithmic derivative of the Hankel determinant 
with respect to t , and (3 n , r n , which is an essential step in establishing the Painleve 
equation. Whilst the difference relations (I4.14p - fl4.16p and f )4.19p look rather complicated, 
these can be manipulated to give us insight into the recurrence coefficients a n and /3 n . To 
this end, summing (I4.14p and (14.151) gives us a simple expression for the recurrence coefficient 
a n in terms of R n : 

a n = 2n + l + a + \-tR n . (4.20) 

From (I4.16P and f 14 . 1 9 1) we find after a minor re- arrangement 

p n {Rn + Rn-i) = fi n - n(n + a) - r n (a + A + 2n). (4.21) 

Now substituting ( H~T9"j) and ( H~2"Tj) into either ([437]) or ( H~T8"j) to eliminate R n and f? n _i 
leaves us the following very simple form for 5Zj=o Rj-> w hich will play a crucial role later, 

n-l 

tJ2 R J =n{n + a + \) - fi n -tr n . (4.22) 
3=0 

But in view of (I4.20p . we also have an alternative representation of J^J=o -^i' name ly> 

n— 1 n— 1 

t^Rj = n(n + a + A) - = n(n + a + A) + p 1 (n). (4.23) 

j=0 j=0 



In summary, we have obtained two different ways to express Y2j R31 an< ^- comparing (14. 2z 
with (I4.23P gives us the important relation (I4.13p . linking Pi(n) to f3 n and r n . We are 
now in the position to find an expression for f3 n in terms of r n and R n . This is found by 
eliminating R n -\ from (I4.2ip and (I4.19P resulting in (14. lip . End of Proof. 

4.2.2 t Evolution and Painleve V 

In the next stage of the development, we vary t and n. The differential-difference relations 
generated here when combined with the difference relations obtained previously will give us 
the desired Painleve equation. A straightforward computation shows that 

j t \ogh n = R n . (4.24) 



But, from fl2TT5]) . it follows that 



(]l - Pn(Rn-Rn-l) (4.25) 

= PnRn- r2n ~ Xr '\ (4.26) 
tin. 
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where the last equality follows from ( 14.19]) . 
Differentiating 

/>oo 

0= / x a (x + t) x e- x P n (x)P n . 1 (x)dx 
Jo 

with respect to t produces 

= A / x a (x + t) x - 1 e' x P n (x)P n ^ 1 (x)dx + 
Jo Jo 

P n _i(x)P n (x) 



x a {x + t) A e 



A„— x 



, / w dLaig (x,t)dx + h n _i —p^n), 



resulting in 



Upon noting (I2.13p . this implies 



d 

dM H) 



da r 



dt 



Now differentiating (14.131) with respect to t and noting (I4.27p . we find 



d 



PiO) 



df3 n d 



(tr r , 



dt dt dt 

The above result combined with (14.261) gives 



dfin 
dt 



r n - 1- 



dr n 
dt 



dr ri 



dt 



dt 



P n -i(x)dx 



(4.27) 



(4.28) 



(4.29) 



We now come to a key Lemma which gives the first order derivative of r n (t) and R n (t) 
with respect to t, and where n appears as a parameter. 

Lemma 5 The auxiliary variables r n and R n satisfy the following coupled Riccatti equa- 
tions, 



dr n r~ - Ar r , 



dt 



Rr, 



1 — Rn 



r n (2n + a + A) + 



Rn 



+ n(n + a) 



and 



Furthermore, 



dR 

2r n = t— - + X-R n (t + 2n + a + X-t R n ). 
dt 



y(t) = y(t,n) := 1 



(4.30) 



(4.31) 
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satisfies the following second order non-linear ode, 



» ( V _ £ 4. ^~ 1 ) 2 _ ^ , (2n + 1 + a + A) y _ y(y + 1) 

* 2y(y - 1) l2/ J t + t 2 \2 V 2y) + t 2(y - 1) ' 

(4.32) 



which is recognized to be a 



' c? a 2 

Pv ( y,-y,2n + l + a + A,-l/2 



Proof Because (14. lip expresses (3 n as a quadratic in r n , we see that r n satisfies the 
Riccatti equation (I4.30p . Eliminating r n+1 from (14.151) and (14.281) . and upon referring to 
(14. lip , we obtain (14.3 ip . Next, we simply substitute r n (t) from (I4.3ip into (I4.30p . to see 
that Rnit) satisfies a second order non-linear ode in t, in which n, a, and A appear as 
parameters. A further linear fractional change of variable 

Rn(t) = 1 - - 1 or y(t) = l 



l-y(t) yw I -Rnit)' 

establishes that y(t) satisfies the Painleve V displayed in the Lemma. End of Proof. 

We begin here a series of computations which ultimately give rise two integral represen- 
tations of the Hankel determinant of interest, i.e., 



D N (t, A) = det x j+k - 2 (x + t) x x a e- x dx^j 



l<j,k<N 

given in Theorems 1 and 2. To this end, an easy computation shows that 

N-l N-l 



H N (t) := t^logD N (t,\) = t^J2 lo Z h i = t J2 R j 

3=0 3=0 
= N(N + a + A) - j3 N -tr N (4.33) 

= N(N + a + X) + p l (N), (4.34) 

where the last two equations follow from (14.121) and (I4.13P of Lemma 4. Integrating (14.331) 
with respect to t , while noting (14.111) . (14.311) and -Rtv(^) = 1 — 1/(1 — y(t)), we obtain the 
result stated in Theorem 1. 

To obtain the second integral representation for D n (t,X) stated in Theorem 2 (i.e., in 
terms of H n (t) ), we note that from (14.271) . (14.331) . and (I4.34p . we obtain expressions for [5 N 
and r N in terms of H N and H' N , 

N = N(N + a + X)+tH' N -H N (4.35) 
r N = -H' N . (4.36) 
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What we need to do is to eliminate to find a functional equation satisfied by H N , H' N 
and H'^. For this purpose, we examine two quadratic equations satisfied by , one of 
which is simply a rearrangement of (14. lip and reads 



N 



Xr N 



R 



N 



+ f3 N R N = p N - r N (2N + a + A) - iV(iV + a). 



(4.37) 



The other follows from a derivative of the first equation of (I4.35p with respect to t and 



PnRn — 



r 2 N - Xr N 



R 



N 



tH'lf. 



(4.38) 



Solving for R N and 1/Rn from the linear system (14.371) and (14.381) . we find 

tH" N - (2N + A + X)r N - N(N + a) 



2R N 
2 

Rn 



1 + tH' N -H N + N(N + a + \) 
-tH N + (t + 2N + a + X)H' N -H N + N\ 



(4.39) 
(4.40) 



N 



where we have replaced (3n and in terms of if at, if^and H N with ( 14. 35ft . The product 
(14.391) and (I4.40H gives us the desired a— form ( 14. 5ft . 

It is finally worth noting that with D]y(t, A) =: t SN D]y, we find, after a little computation 
that D N satisfies the Toda molecule equation [21] 



(4.41) 



d 2 ~ Dn + iDn-i 

-—logV N = ~ . 

dt 2 & D 2 N 



4.3 The Discrete a— Form 

As an alternative to the Jimbo-Miwa-Okamoto a— form satisfied by in Theorem 2, 
we state (without proof) a non-linear difference equation satisfied by H N . We believe this 
to be new. Specifically, we find that the logarithmic derivative of the Hankel determinant 
in (13.161) . generated by the deformed Laguerre weight in (13. 14j) . satisfies a second order 
non-linear difference equation, which we call the discrete a— form, 



iV(iV + a)t + (A 2 H N + t)(H N - 5 



-i 2 



A' J 



A 2 H N + 2N + a + A + t 



A 



N(N + a)t + (A 2 H N + t)(H N - 5 



N) 



N 



H N + 



N(N + a)t + (A 2 # 



N 



A 2 H N + 2N + a + A + t 



A 2 H N + 2N + a + \ + t 



N+l 



Hn){Hn — H 



N-l, 



where A 2 H 



N 



Hn+i — H 



N-l 



The initial conditions are Hi(t) = -^logDi(t, A), H 2 (t) 



dt 



\ogD 2 (t, A) with Di(t,X) = Ho(t), D 2 (t,X) = j^o(t) ^(t) — /J, 2 (t), and the moments are 



defined in (0071) . 
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We do not further elaborate on this result, but present it here since we believe that it may 
provide a useful tool for efficiently computing in an iterative manner, and give further 
insight into the mutual information of a single-user MIMO system. A thorough analysis of 
the implications of this non-linear difference equation in regards to wireless communication 
applications is the subject of on-going work. 

5 Concluding Remarks 

This is a short review is about Hankel determinants that arise in the information-theoretic 
study of MIMO communication systems. We consider the capacity of single-user MIMO sys- 
tems, in which case the determinant of interest is generated from a certain deformed Laguerre 
weight. We have obtained two exact integral representations for this Hankel determinant, 
one of which is described in terms of the a— form of a particular Painleve V differential equa- 
tion. We have also stated an alternative representation, involving a second-order non-linear 
difference equation. 
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